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Abstract. We introduce a class of orthogonal polynomials in two variables which generalizes 
the disc polynomials and the 2 -D Hermite polynomials. We identify certain interesting 
members of this class including a one variable generalization of the 2 -D Hermite polynomials 
and a two variable extension of the Zernike or disc polynomials. We also give ^-analogues 
of all these extensions. In each case in addition to generating functions and three term 
recursions we provide raising and lowering operators and show that the polynomials are 
eigenfunctions of second-order partial differential or g-difference operators. 
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1 Introduction 

One of the earliest orthogonal polynomials are the Legendre polynomials which are orthogonal 
with respect to dx = (linear measure) on the unit interval. They are the spherical harmonics 
in M 2 . The ultraspherical polynomials are the spherical harmonics on M m and are orthogonal 
with respect to (-*- — x 2 ) u dx on the unit interval. If we replace x by x/y/v and let 

v —> oo the measure becomes 7r _1 / 2 e _x " i dx and the ultraspherical polynomials, properly scaled, 
tend to Hermite polynomials. The next sequence of polynomials in this hierarchy are the Jacobi 
polynomials which are orthogonal on the unit interval with respect to 2 r(a+i)T(/3+io~ 2 ^ — 
x) a (l + x)Pdx. The Laguerre polynomials arise when x is replaced by —1 + 2 x/a and a —> oo. 

The 2D history is somewhat parallel. In the late 1920’s Frits Zernike was working on optical 
problems involving telescopes and microscopes. He introduced polynomials orthogonal on \z\ < 1 
with respect to the area measure, so this is the 2D analogue of Legendre polynomials, see [38, 39]. 
The more general disc or Zernike polynomials are orthogonal with respect to (1 — x 2 — y 2 ) v dxdy 
on the unit disc. Again if we replace (x,y) by (x/y/u, y/y/u) and let v —> oo this measure, 
properly normalized becomes e~ x ~ y dxdy on M 2 . The polynomials become the 2ZJ-Hermite 
polynomials introduced by Ito in [22] in a different way. They are defined by 

mAn / \ / \ 

H m ,„(z 1 ,z 2 ) = J2 (”) (")(-i)‘*k"-‘(f)"-\ 

*This paper is a contribution to the Special Issue on Orthogonal Polynomials, Special Functions and Applica¬ 
tions. The full collection is available at http://www.emis.de/journals/SIGMA/OPSFA2015.html 
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where m An = min{m, n}. In 1966 D.R. Myrick [27] considered only the radial part of the 
polynomials orthogonal with respect to the Jacobi type measure (x 2 + y 2 ) a ( 1 — x 2 — y 2 Ydxdy on 
the unit disc, see also [36]. Addition theorem for the disc polynomials were found by Sapiro [29] 
and Koornwinder [24]. Maldonado [26] very briefly considered a class of radial functions which 
can be used to generate 2D orthogonal polynomials. P. Floris [10] introduced a (/-analogue of 
the disc polynomials proved an addition theorem for them, see also Floris and Koelink [11]. 

There are several other families of orthogonal polynomials in two and higher variables. Koorn- 
winder’s survey article [25] covers most of what was known up to the time of its publication in 
1975. In the 19th century Hermite found families of orthogonal and biorthogonal polynomials, 
see [14, 15], [8, Chapter 12] and [3, 4], The more recent book [7] contains a treatment of the 
theory and applications of orthogonal polynomials in 2 and higher variables. One recent related 
reference is [33]. The 2U-Hermite polynomials have may applications to physical problems, see 
[2, 6, 30, 32, 34, 35]. Mathematical properties of these polynomials have been developed in 
[12, 13, 16]. A multilinear generating function, of Kibble-Slepian type, is proved in [18]. In [21] 
we gave two (/-analogues of the 2D-Hermite polynomials and studied the polynomials in great 
detail. We also introduced a different (/-analogue of the disc polynomials. 

In this paper we first identify a general class of two variable polynomials whose measure is 
the product of the uniform measure on the circle times a radial measure. This class not only 
include the 2U-Hermite polynomials and their (/-analogues but it also contains the Zernike (or 
disc) polynomials and its (/-analogues. This will be done in Section 2. One special feature of 
this class is that multiplication by z or z results in a three term recurrence relation. This special 
feature is unlike what happens for general real 2 — D systems [7]. In fact the earlier work [18] 
motivated Yuan Xu [37] to reformulate the theory of orthogonal polynomials in several variables 
in terms of complex variables. 

In Section 3 we give a one parameter extension of the 2D-Hermite polynomials. They are 
2 D analogues of Laguerre polynomials and will be denoted by {Zm,n{z,z)}. We record their 
definition, orthogonality relation, and the three term recurrence relations in Section 3. In 
Section 3 we also derive several differential properties of these polynomials. It must be noted 
that the combinatorics of these polynomials have been explored in [19]. Section 4 contains 
a treatment of a two parameter family of disc polynomials. They are orthogonal with respect 
to {x 2 + y 2 ) a ( 1 — x 2 — y 2 )^dxdy on the unit disc. They were briefly studied in [20] but here we 
derive new generating functions and study their differential-difference properties and show that 
they satisfy a Sturm-Liouville system. 

The g-analogues start with Section 5 where we introduce two ^-analogues of the 2D-Laguerre 
polynomials {Z^n{z, z\q)}. Both can be considered as a (/-extension of the 2D-Laguerre polyno¬ 
mials, one for q £ (0,1) and one for q > 1 . As expected the case q > 1 leads to an indeterminate 
moment problem while the case q £ (0,1) is a 2D-extension of the Wall or little (/-Laguerre 
polynomials, [5, 23]. A (/-analogue of the two parameter 2D-disc polynomials is introduced in 
Section 6 and many of their properties are developed. In Section 7 we find polynomial solu¬ 
tions of a second-order partial differential equation. Section 8 contains a system of biorthogonal 
polynomials which is a possible 2U-analogue of the Askey-Wilson polynomials. 

We shall use the notations 


(D q , z f)(z) 


m - f(qz) 

(i -q)z 


(<W)(*) 


a , 


1 - q a 

1 ~q' 


mV n = maxjm, n}, 


z(D qtZ f)(z), 
m An = min{m, n} 


and 



= zd z . 


d : 


5 ; 
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Throughout the paper we will use Z\ and Z 2 - The polynomials we consider here are polyno¬ 
mials of two real variables, so we shall always assume that Z 2 = z\. 


2 General constructions 


We now give a concrete construction for general 2D systems. Given JVgNU { 00 }, let v{0\ N ) 
be a probability measure on the circle T = {e ie | 0 < 6 < 27t} of Fourier type 

[ e^ m -^ e dv{ 6 - N) = S m , n , m,n = 0,l,...,N 
J T 

and { 4 > n (r-,a)} be a system of orthogonal polynomials satisfying the orthogonality relation 

POO 

/ 4 > m (r]a)(j) n (r;a)r a dn(r) = ( n (a)5 m ,n, a> 0. 

Jo 

It is assumed the positive Borel measure /j. does not depend on a. Let 


4> n {r-a) = ^ ~2 c j(n,a)r n j , Cj(n,a ) £ 
3 =0 


and define polynomials 

^2) P) — 

This defines the polynomials for all m, n and z\ . Z 2 £ C. It is clear that 


z™ n <t>n(z\Z 2 ',m — n + ft), m>n, 
fn,m(z 2 ,Zl\ /3), m < 71 . 


fm,n{z,z;P) = fn,m(z,z]/3), m < n. 
From (2.2) it is clear that for m > n 

Z\fm,n(z\i ^2) ft T 1) — ./m+1 ,n(zi , Z2-, ft ), 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


and 


d 

~'i-QQfm, n (z, Z\ ft) 


{m-n)f m , n (z,z;ft). 


Therefore, 


{$z i - Sz 2 )f m , n (zi,Z 2 ]ft) = {m- n)f m>n (zi,Z 2 ift), m > n. 

Similarly, 

(Sq, Zl - q m ~ n 8 q , Z 2 )fm,n{zi,Z 2 -,ft) = [m - n\ q f m>n (z 1 ,z 2 ;ft), m > n. 

It must be noted that f- m ,n{z \, 22 ) is defined for z\,Z 2 £ C but the next theorem gives an 
orthogonality relation on M 2 . This is exactly analogous to the classical univariate orthogonal 
polynomials which are defined for z £ C but their orthogonality is on a subset of M. 

Theorem 2.1. Given N £ N, for nonnegative integers m, n, s, t such that m+t < N, n+s < N 
the polynomials {f m ,n(zi, 22 ; ft)} satisfy the orthogonality relation 

/ fm,n(z,z\ft)f S f{z,z\ft)dv{Q\N)dii{r‘ 2 -,ft) = CmAn(\m - n\ + ft)S mjS 6 n>t , (2.4) 

Jr 2 

where dp,(r; ft) = r^d/i(r). If N = 00 , (2.4) holds for all nonnegative integers m, n, s, t. 
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Proof. There is no loss of generality in assuming m > n. Otherwise, we simply apply (2.3) 
to switch the orders of m and n, s and t , then consider the consider the complex conjugate of 
the obtained integral instead. By performing the integral over 9 we see that the left-hand side 
of (2.4) equals 

/ fm,ni z ,z',P)fs,tiz,z-,/3)di>i6)diJ,ir 2 ;l3) =0 
Jr 2 

unless m + t = n + s. Let a = m — n = s — t. Therefore the above integral equals 

roc n t 

/ Y, Cj (n, a + p) r 2n + a ~ 2 i Y c kit, « + (J)r 2t+a ~ 2k r 2 ^dfi(r 2 ) 
j=0 k=0 

/*00 n t 

J ^2 c,-(n, a + P)r n ~ J Y Cfc ^’ a + /3)r*“ fc r a+/3 dn( r ) = Cn(a + P)6 n ,t, 


3=0 

and the proof is complete. 


k =0 


Remark 2.2. For any positive integer N, if we take N ) as the discrete probability measure 
used in FFT (fast fourier transform), we get finitely many 2D orthogonal polynomials. But in 
most of our applications we use N = oo and the circle measure ^ for dv{6) unless specifically 
stated. 

We now come to the three term recurrence relation. 

Theorem 2.3. The polynomials {f m ,n{z i, ^ 2 ; /?)} satisfy the three term recurrence relation 
r , o\ coin, m-n+l + P) 

Z2fm+\,n\Zx,Z2\P) = „ ^ , , - , fm+l,n+l{Zl, Z 2 ] P) 


co(n + 1 , m — n + ff) 
c 0 (n,m — n + 1 + /3)c n+ i(n + l,m - n + /3) 
co{n + 1 , m — n + /3)c n (n, m — n + /3) 


fm,n ("i 1 ; ft) > m>n. 


(2.5) 


Proof. Since the polynomials {^) n (r; a)} are orthogonal on [0,oo), their zeros are in (0, 00 ). 
Hence ^> n (0; a) 7 ^ 0 for all n and all a > 0. From the Christoffel formula [17, Theorem 2.7.1] we 
see that 

rcfnir;a + l) = ^ —r[<?Wi(D a)</>n(0; a) - <£ n (r; a)<£ n+ i(0; a)], 


cq in + 1 , a)(f> n i0; a) 


that is 


rMn°+0= C °i"’ a , + 1 ! fe+,(r ;a )- C ° ( "-r + 1 , )C "t l( r + 1 , ,a> ^(r; a ). ( 2 . 6 ) 


coin + 1 , a) 


co(n + l,a)cn{n,a) 


Now take a = m — n + (3, replace r by Z 1 Z 2 then multiply both sides by z\ m n to estab¬ 
lish (2.5). ■ 

One can also prove that [20] 

Z\ fm.ni%l, Z 2 , /3 s ) t ' ~ : wr fm-\-l,nizi > z 2 1 P) — 1 (^1 , z 2', P '), (2.7) 


where 




coin, m + 1 — n + (3) ‘ 

coin, a + l)ci(n, a) — coin , a)ci(n, 1 + a) 
co(n — 1,1 + a)co(n, 1 + a) 


and a = m — n + (3. 
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It must be noted that the 2D-polynomials may have recurrence relations containing more than 
three terms [7]. It is surprising that the polynomials in the special class we are considering have 
three term recursions. Another important point is that the definition of the polynomials f m ,n 
indicates that it has either a trivial zero at 2 = 0 or the zeros will lie on circles whose radii are 
the square roots of the zeros of ^> n (r; m — n + /3). The zeros of <^> n (r; m — n + /3) are positive, 
real and simple. 

Let 


rcj) n (r] a) = a n (a)c/) n+1 (r; a) + c n (a)(/) n (r; a ) + b n (a)</> n -i(r; a ) (2.8) 

be the three term recurrence relation satisfied by the polynomials {( f) n (r;a )}, for some real 
numbers a n (a), b n (a), c n (a). By matching the corresponding coefficients we find that 

coin, a) 


®n(®) — 

b n (a) = 

Cn(a) = 


c 0 (n + 1, a) ’ 

co(n, a)c 2 (n, a) — ci(n, a) 2 coin, a)c 2 (n + 1, a) — c\{n, a)ci(n + 1, a) 


coin — 1, a)co(n, a) 
ci(n, a) ci(n + 1, a) 


co{n — 1, a)co(n + 1, a) 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


coin, a) c 0 (n + l,a)‘ 

Making use of (2.1) and (2.2) we conclude that the recursion (2.8) becomes 
[ziz 2 - c n (m - n + /3)\f m ,n(zi, Z 2 \ f3) 

= a n (m - n + /3)f m +i,n+iizi, Z 2 ', /3) + b n im -n + P)f m -i,n-iizi,Z 2 -,P). (2.12) 

Example 2.4. Consider the case of the 2 -D Hermite polynomials when the coefficients in (2.2) 
are given by 

n!(n + a)!(-iy 

Cj{n,a) = —7-v-77--T7. 

j!(n - j)!(n + a - j)\ 

The recursion (2.5) leads to the second three term recurrence relation 


Z\H m n iz\, Z2) — ^H m n —\iz\, Z2 ) T TLm+l.n(^lj Z2), 
Z 2 H m , n (zi,Z 2 \q) = mH m -i t nizi, Z 2 ) + H m . n+ iizi,Z 2 ), 


while the first can be proved by direct computation. 

Similarly the case when the coefficients in (2.2) are given by 


Cjin , a) 


iq',q)niq-,q)n+ai- i) J 
iq',q)j{q\q)n-jiq\q)n+ a -j 1 


leads to the recurrence relations for this ^-analogue, see [21]. 


Remark 2.5. If we demand that both (2.2) and (2.3) are valid for all m, n, z, then we must 
have 


Cjin, a) 


fn\ a jin, a) 
\j) (n + a-j)\ 


for the classical cases, and 


Cjin, a) 


n\ a jin, a\q) 
j J n ( q'i q)n+a-j 


(2.13) 


for the (/-analogues, where Cjin,a) are the coefficients in (2.2). These factorizations imply that 
</> n (r; a) must be a generalization of Laguerre orthogonal polynomials. However, in most of the 
cases, (2.2) is only valid for m > n, whereas (2.3) acts as a kind of analytical continuation to 
the case m < n. 
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Remark 2.6. In view of the defining equation (2.2) the question of finding the large rn, n 
asymptotics of the two variate polynomials f m ,n is equivalent to finding the large m, n behavior 
of i j) n (r,m — n + j3). It will be of interest to carry out this program at least for some special 
systems, including the Freud type polynomials orthogonal with respect to x a exp (-p(x)), where p 
is a polynomial with positive leading term. 

We now consider an inverse to (2.6). 

Proposition 2.7. Given a nonnegative measure on (0,oo) 

d/i(x; a) = x a dv{x), x > 0, a > — 1 

such that 

x n dg>(x] a) < oo, n = 0 , 1, .... 

Let {cj) n (x; a)}™ =0 satisfy 

fm(x ) Ol)(j)n(x , tt)(i/x(x, tt) — Cn(®)^?n,n- 

Then 




4> n (x- a) - a n (a)(f> n (x; a + 1) = b n (a)(j) n - i(x; a + 1), 


(2.14) 


where 


b n {a ) = 

Ojyiiof) = 
Furthermore, 


n (x-,a) = y jCj (n,a)x n - j , 
j =o 

co(n, a + l)ci(n, a) — co(n, o)ci(n, a + 1) 


co(n — l,o + l)co(n, a + 1) 
c 0 (n,o) 


c 0 (?z, a + 1)' 


POO 

/ x<f>n(x] a)(t> n - i(x; a + l)d/r(x; a) = b n (a)( n -i(a + 1). 

Jo 


(2.15) 


Proof. It is clear that cj) n (x-,a) — a n (a)(j) n (x-, a + 1) is a polynomial of degree at most n — 1, 
then 


n—1 


(j) n (x-, a) - a n {a)(t) n {x\ a + 1) = ^ ajf>j(x; a + 1), 

i=o 


where 


a 


POO 

<(j(a + 1) = / {4> n (x-a) - a n (a)4> n (x;a + l)}4>j(x; a + l)d^{x\a + 1) 

Jo 


= / f>n(x-, a)(fj(x] a + l)d/u(x; a + 1) 

Jo 

P OO 

= / xcj) n (x; a)(fj(x; a + l)d^(x; a). 

Jo 

Clearly aj = 0 for 0 < j < n — 2 and o n -i = b n (a). Then b n (a) is obtained by matching the 
leading coefficients of two sides. ■ 
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Corollary 2.8. Under the assumptions of Proposition 2.7, let 


4 > n {x\ a + 1) = ^ A j(n, a)4>j(x ; a). 
j =o 


Then 


\ ( \ ( —J TT b n -k(a) . 

Xj(n,a) = \ J_| -- 777, 0 < j < n — 1, 


and 


A n {n,a) = 


Cn(a) 


a j { a ) a n-k{c*) 

1 


Q>n (c *0 


n— 1 

n 


b n -j(a + j)c 0 (n - j, a + j ) 


Co(« + n) a n -j-i(a + j)c 0 (n -j-l,a + j)' 


Furthermore, we have 


71—1 


d n (a)(j) n (x;a) = <f> n (x\a + 1) - a; ^ fj(n, a)4>j{x-, a + 1) 

3=0 


and 


n —1 


<l>n(x;a) = ^> n (0; a) - x Y a)</>j(x-, a + 1), 

3=0 


where 


d n (a) = 


Cn(n, a + 1) 
c n (n,a) 


fk{n,a) = Y 


3=0 


Aj(n, a)Xj(k, a)(j(a + 1) 

Cfc(a + 1) : 


0 < < n — 1, 


/ X (-l) fc Co(a)co(0, a + l)c n (n, a) k j4 b k _j(a ) 
g k (n,a) = - A ^ ^ -_[]_ 


Cfc( a +1) 

Proof. From (2.16) and (2.14) we get 


3=0 


® k—j (®) 


n—1 


i i> n (x\ a) — a n (a) Y^ Xj(n, a)<j>j(x; a) = b n (a ) Y^ Xj(n — 1, a)4>j{x\ a). 

3=0 3=0 


Then 


1 - A n (n, a)a n (a) = 0, 


and 


(2.16) 


(2.17) 

(2.18) 


(2.19) 


—a n (a)Xj(n, a) = b n (a)Xj(n - 1, a) 
for 0 < j < n — 1, which lead to (2.18) and (2.17) respectively. 
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From (2.15) we get (2.16) we get 

POO 

b n (ai)(n-i{at + 1) = A n _ i(n - l,a) / xc/) n (x; a)^„_i(x; a)dn(x\ a) 

Jo 

x(j) n (x ; a)i?i ri _i(x; a)dfi(x ; a). 


1 


®n—l(o0 io 

Since {</> n (x; a)}^T 0 are orthogonal polynomials with respective to the measure dfi(x-,a), hence 
X<f>n—1 i^'i oJ) = (s?! Ql) T 5 n (f> n —l(x] <a) T C n (f>n —2 (®j Oi'), 

where 

coin — 1, a) 


Vn = 


c 0 {n, a) 


Then 


M a )Cn-i(a + 1 ) = 


VnCn(cv) 

a. n _i (q) 


and 


_ b n (oi) co{n, a) ^ ^ ^ 

Cni a ) — r \ ( 1 \ Cn—l (c^ T 1) 

a„_i(a) coin - 1, a) 


n— 1 


K-j (a + j ) c 0 (n - j, a + j ) 


Co(a + n ) n o^-j-r(a + j) c 0 (n - j - 1, a + j ) : 


which gives (2.19). 
Let 


n—1 


4>nix', a + 1) - d n ia)4>nix- a) = x ^ /j(n, a)^(x; a + 1). 

i=o 

For 0 < /c < n — 1 we have 

roc 

/ {0n(x, a + 1) - d n (a)^n(x, a)} </> fc (x, a + l)d/r(x; a) 

Jo 


a)4>jix, a + 1) 



> (j)kix, a + l)d/r(x; a) 


= fkin,a)( k ia + 1), 


which is 


r oo 


fkin, a)Cfc(a + 1) = / </) n (x, a + l)^ fc (.x, a + l)d/i(x; a) 

Jo 

k 


'^2\ j in,a)\jik,a)(jia + 1 ), 
3=0 
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or 


A j(n, a)Xj(k , a)(j(a + 1 ) 


Cfc( a + 1) 


f k {n,a) = Y.-^’'' 3yn ^ 

3 = 0 

On the other hand, let 


71—1 


(/>n(x;a) = cj> n (0] a) - x^2gj(n,a)(f)j(x;a + 1), 

3=0 

then, for 0 < A: < n — 1, we have 

POO 

/ {4>n(x; a) - </>n( 0 ; a)} <f> k (x, a + 1 )dg(x; a) 

Jo 

n -1 'j 

® gj(n, a)4>j(x-, a + 1 ) > 4> k (x, a + 1 )dg(x; a ) 


j=o 

n—1 


y, gj{n , a) 0 j(x; a + 1 ) > 4> k (x,a + l)dg(x;a + 1 ) 

3=0 


= - 9 k{n,oc)C k (a + 1 ), 


that is 


POO 

g k (n,a)( k (a + 1) = / 0 n (O; a)cf k {x, a + l)dg(x; a) 

Jo 


i o 

= 4> n ( 0; a)Xo(k, a)<j> o(0, a)Co(a) = c n (n, a)c 0 (0, a)Co(a)Ao(A:, a) 

fc_1 b k -j(a) 


(-l) fc Co(a)co(0, a + l)c n (n, a) J j 


/=0 a k-j(ot) 


that is 


, \ (-l) fc Co(a)c 0 (0,a +l)c n (n,a) -r-r 6 fc -j(a) 

g k {n,at) = -- 7 ——rr- M - 7 -^. ■ 

Cfc(a + 1) 

Corollary 2.9. Let 4> n (x, a), a n (a), b n (a ) as in Proposition 2.7, form > n— 1 file equation (2.7) 
is in file /orm 

2l/m,n(zi,Z2) ~ Vm,n /m+l,n(^l ? ^ 2 ) — 'U j m,nfm,n—1 (zi,z 2 ), (2.20) 


where 


Um,n = b n (m-n), v m>n = a n {m - n). 


Proof. For m > n — 1, applying the definition f m ,n(zi, z 2 ) = z™ n (j) n (ziz 2 -,m — n), m > n to 
get (2.20). ■ 


We next consider differential or ^-difference equations satisfied by f m ,n(zi, z 2 ). 
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Theorem 2.10. Assume that cj) n (r;a) satisfies the second-order differential equation 
A n (r, a)6 2 f + B n (r, a)5 r f + C n {r , a)f = 0, 

where 

$rf = r-jrf- 

or 

Then f m ,n(zi ■< z 2 ; fi), for a = m — n + fi and rn > n, satisfies the second-order partial differential 
equations 

A n (z 1 z 2 , a)5 2 Z2 f + B n (ziz 2 , a)5 Z2 f + C n (z 1 z 2 , a)f = 0, 

A n (z 1 z 2 ,a)8 2 J + {B n (ziz 2 ,a ) - 2aA n (ziz 2 ,a)}5 zl f 

+ {a 2 A n (ziz 2 ,a) - aB n (z 1 z 2 ,a ) + C n (ziz 2 ,a)} f = 0. 

Similarly, assume that </> n (r;a) satisfies the second order q-difference equation 

A n (r, a)8 qr f + B n (r, a)5 q ^ r f + C n (r , a)f = 0, 

where 

8q.rf = fD qr f. 

Then f m ,n(zi, z 2 ), for m>n and a = m — n + fi , satisfies the q-partial difference equations 

A n (z!Z 2 , a)5 qZ2 f + B n (ziz 2 , a)5 qtZ2 f + C n (ziz 2 , a)f = 0, 

A n (ziz 2 ,a)8 2 zi f + { q a B n (ziz 2 ,a) - 2[a\ q A n (ziz 2 , a)}5 q , Zl f 

+ {[a] q A n (ziz 2 , a) - \a\ q q a B n (ziz 2 ,a) + q 2a C n (ziz 2 ,a)} f = 0. 

Proof. Observe that 

<? > >n(ziZ2 ) = Z^ a f m ,n (zi,Z 2 ), 

$ziZ2 ( t > n(z\Z2) — Zi fi Z 2fm,n(zi,Z2') — Z^ {^z± f 

^zfZ 2 ^n[zxZ 2 ) = Zi a 8 2 Z2 f m , n (zi, z 2 ) = z^ a {6^ - 2 aS Zl + a 2 } f 

fiq,ziZ2 ( t > n(z\Z2) — Z ^ $q,Z2 /m,n(^l > Z 2 ) — Q Z^ {$q,zi 

6 lz lZ2 ^n{z 1 Z 2 ) = zf a 5 2 qtZ2 f m , n (z 1 ,z 2 ) = q~ 2a zf a [5 2 Zl - 2[a] q 5 q , zi + [a} 2 }f. 

Then the proof of theorem follows from substituting the required combinations into the corre¬ 
sponding differential equation or ^-difference equation for cf>(r,a). ■ 

A simple application of Fubini’s theorem establishes the following theorem concerning gen¬ 
erating functions. 

Theorem 2.11. For z\, z 2 ,u,v € C such that 


2>2 ~T~ ls r j(Ai22,ra)| < °o, 


1=0 


E 

i=o 


\UZlY 

(?;?): 


\g j (z 1 z 2 ,uv)\ < oo, 


where 


g Q (x, z) = ^2 Tk(x; a)z k . 


k =0 
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Then 


E 

m>ri> 0 

E 


(m — n)\ 


( Q> Q)m—r 


fm,n( z li Z 2 ) = Y ^~~r 9j( z l z 2,Uv), 


1=0 


m>n> 0 

Similarly, for z\,z 2 ,u,v £ C such that 


fm,n{zi, Z 2 \q) = Y 7~~4 9j{ziz 2 ,uv). 
j=0 ' q ' q >3 


Y\ uz i\ j \9j(ziz 2 ,uv)\ < 00 , 
i=0 


then 

OO 

Y U m V n fm,n(z\,Z2 ) = ^(uz+ffj (ziZ 2 , Uv ), 
m>n>0 j=0 

00 

^ M m r> n / m , n (zi,z 2 k) = 

m>n>0 j=0 

Theorem 2.12. Assume that 
w a+/3 (x) = x a wp(x) 


( 2 . 21 ) 


and 


+ + a) = — —d™ (w a (x)x n ), n = 0,1, 
w a (x) 

Then for m> n we have 

t „ . a ^_ f Zi\ m ~ n df 2 (' w p ( Zl z 2 )z?) 

fm,n\ z l,Z 2 ,f3) — ( ) / \ j 

\Z 2 J Wp{ZlZ2) 


d. 


^2 


/m+l,n+l(^-1) ^-2! /3) — 
fm,n(z\, Z 2 \ /3) — 


m+l—n 


W / 3(z 1 Z2)fm+l,n(zi,Z 2 -, f3) 


W 


'p(ziz 2 ) 


d Zl {wp(ziZ 2 )z 


fm+l,n+l(zi, Z 2 ] /?) — 


Wp{ziz 2 ) 

8 Z1 (Wi3(ziZ2)f m +l,n(zi,Z 2 ]/3)) 

w / 3 (z 1 z 2 ) 


and 


Z2 


m+l—n 


d Z2 ( ( — ) Wp(ziZ2)fm+l,n(zi,Z2\P) ) = 9 Zl (wp(z 1 Z 2 ) • 

. %1, 


( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 


Furthermore, if for some nonnegative integers s, t such that 

w' (x) ^ ajxj 

w /3\ x ) 3=0 

1=0 
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where {ajR =0 and {6j}* =0 are certain real numbers. Then form m > n we have 


j =o 


i d Z2 {wp{z 1 Z 2 )f m +j-\, n { z ^ z 2 \ P)) 


3~2 


W 


] p{ziz 2 ) 


'p ] a j z 2fm+j,n{zi, Z 2 ',P), 
3=0 


< d 

i2 b r 

3=0 


21 V "1 


n—m—j -\-1 


'Wp(ziZ2)fm+j-l,n(zi,Z2; 


n—m—j 1—1 / \ 

z 1 J z 2 J wp(z 1 z 2 ) 


^ ^ a j Z 2 fm+j,n { Z 1 1 z 2 1 P) 
3=0 


and 


3=0 


,. z 3 

3-2 


d Z2 {wp{ziZ 2 )f m +j-l,n{.Z\,Z 2 ] f5)) 
wp(ziz 2 ) 


* d Zl (Zi m j+1 Wp(ziZ 2 )f m +j-l, n {zi,Z2-, 

E l >: 


n—m—j I—? / \ 

z 1 J Z 2 J wp(ziz 2 ) 


3=0 


Proof. From 


d n Z2 [w a (ziz 2 )(ziz 2 ) n ] = dl (wpiz^Xz^D = (wp(z lZ2 )z^) 

= z?w a (z 1 z 2 )(/) n (z 1 z 2 ; a) = z'fzY n wp{ziZ 2 )fm,n{z\,z 2 ]^) 


and 


d 


Z 2 


777.+1 — n 


wp(ziz 2 )f m+ i t n(zi, z 2 ;P) 

( \ m—n 

^) Wp{ziZ 2 )f m +l,n+l{z\,Z 2 \P) 


to get (2.22) and (2.23), (2.24) and (2.25) are proved similarly. Equation (2.26) is 
from (2.23) and (2.25). 

From 


wp(z l z 2 )f m>n (zi,z 2 -, /3) = < 9 ™ (wp(z 1 Z 2 )z^ 1 ) 
we get 

d Z2 (wp(z 1 Z 2 )fm,n{zi,Z 2 ]P)) = d zi (w£\ziZ 2 )Z™ +1 ) . 

On the other hand 

t S 

Y bjz{ +m z{w^\ziz 2 ) = Y a 3 z i +m 4 w p(ziz 2 ) 

3=0 j =0 

implies 

Y b i Z 2 d zi( Z l +mW ?( ZlZ2 ^ = ^2 a 3 Z 2 d z 1 { Z l +mW d( Z 1 " 2 », 

3=0 3 =o 


(2.27) 

(2.28) 

(2.29) 


obtained 


which gives (2.27). Equation (2.28) is proved similarly. (2.29) is obtained from (2.27) 
and (2.28). ■ 
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The following theorem can be proved similarly. 

Theorem 2.13. Assume that 

w a+/3 (x) = w Q+ p(x\q) = x a wp{x\q) 

and 

4>n(x; a ) = (j) n (x ; a\q) = — D ” (w a (x\q)x n ), n = 0 , 1 ,.... 

w a {x) q ’ 

Let z 2 -, p\q) = fm,n(zi,z 2 ;P), then for m > n we have 


z 2 ; /3\q) = 


dq,z 2 {wp{z lZ2 \q)zf) 


d, 


wp(ziz 2 \q) 

m+1 —n 


q,z 2 


/m+l,n+l(^l5 z 2\ — 


wp{z\z 2 \q)f m+ i, n (zi,z 2 -,P\q) 


wp(ziz 2 \q) 


fm,n{zi,Z 2 \P\q) = 


dq , z ! 


/m+l,n+l(^l,^ 2 ;^k) = 


wp(ziz 2 \q) 

+ 21 ( wp(ziz 2 \q)f m +l,n(zi , Z2j /?|g)) 


and 


3, 


q,z2 


-) 

zij 


Z2 


\ m+1 —n 


wp(ziz 2 \q)f m+ i !n (zi, z 2 ; f3\q ) 


zi 


9q,z i + / 3 + ^2|g)/m.+l,n(^l,^2 + |g)) • 


Furthermore, if for some nonnegative integers s, t and real numbers {aj}j =0 
have 


+ ajxi 

D q ,xW0(x\q) = j=0 
3=0 


t/ien for m > n 


+21 + ^ ^(ziZ2|g)/m+j-l,n(+, ^2! /+)) jf ( 

6 i- n—m—j 1-3 -7-n- = E a j4fm+j,n( z U 


3=0 

t 


n—m—j i—7 / i \ 

+ •'w / 3 (ziz 2 |g) 


3=0 


ri, m 7 ^2( w /3( ;gl - S2 k)/m+3- 1 .n(^l,^2;^|g)) _ ■+ j f , m , 

l^ b l Z 2 U+Zi + g) - 2^ a J Z 2/m+3,«(.^l) 2; 2,P|g) 


3=0 


3=0 


and 


h d i’ z m J+lw h( z ^\ ( l)fm+j-i,n{zi,z 2 -,P\q)) 

3=0 


3=0 


'3 z 2 


n—m—j 1—7 / l \ 

: +22 +/+iZ2|g)/m+3-l,n+L, ^2! /%)) 


u+zi+g) 


and {bjYj =0 we 


Z2;P\q), 
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We now discuss monotonicity of zeros of {f m ,n{z, z)}. The following version of Markov’s 
theorem is Theorem 7.1.1 in [17]. 

Theorem 2.14. Let {p n (x;r)} be orthogonal with respect to da(x;r), 
dp,(x;r) = p(x-,T)d/3(x), 

on an interval I = (a, b ) and assume that p(x; r) is positive and has a continuous first derivative 
with respect to r for x E I, r £ T = ( 71 , 72 ). Furthermore, assume that 

x j p r (x;T)dfi(x), j = 0, 1 ,... ,2n- 1 , 

converge uniformly for r in every compact subinterval of T. Then the zeros of the orthogonal 
polynomials p n (x; t) are increasing ( decreasing ) functions ofr, r £ T ; if d{\n p(x - , t)} / dr is an 
increasing ( decreasing ) function of x, x G I. 

The standard version of Markov’s theorem is when /3(x) = x, see [31]. 

By applying Theorem 2.14 to the polynomials {</> n (r; a)} defined in (2.1) it follows that the 
zeros of </> n (r;a) increase with a,a > 0. Now fix n and assume that m > n. In view of the 
definition of the polynomials {f m ,n(z, z)}, m > n, in (2.2) we see that if ?’o is a zero of (j> n (w, a) 
then fm,n{z,z) = 0 for all z on the circle \z\ = y/rfi. Since a = m — n the following theorem 
holds. 

Theorem 2.15. For fixed n the radii of the circles forming the zero sets of fm,n(z,z) increase 
with m for all m > n. 

3 The polynomials {Z^ n (zi, z 2 )} 

Motivated by the class of general 171 systems in Section 2 we define polynomials {Zm}n{zi, £ 2 )} 
by 



Z m]n( z l,Z 2 ) 



(fi + l) m 
(P + 1 )m-k 


(- 1 ) 


n—k „m—k jn—k 


'1 


£0 


for m > n. When m < n the polynomials are defined by 

Z m]n( z i’ 2 2 ) = Z^l(z 2 , zi) 


(3.1) 


Here f3 > — 1. 

These polynomials arise through the choice <j> n {x\ a) = Li? + ^\x), where L[f\x) is a Laguerre 
polynomial [17, 31] 

j (a) / \ __ ( a + l)n fn\ {~x) k _ ( a + l) n ( n \ ( ~x) n ~ k 

n{) n\ 2^\k)(a + l) k n\ ^\k){a + l) n _ k 

Indeed, for m > n 

zW n {zi,z 2 ) = z?- n Li?+ m - n \z lZ2 ) 

and 


ziZ^\zi,z 2 ) = Z^l ljn ( Zl ,z 2 ). 
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It is clear that when /3 = 0 we see that 

H m ,n(zi,Z 2 ) = (-l) n (rn An)\Z$ n (z 1 ,z 2 ). 
Therefore in the notation of Section 2, we have 


Cn( a ) — 


r(a + j3 + n + 1) 


Cj(n, a) = 


(—l)"-J'(a + 0 + l) r 


(3.2) 


nl ’ {n — j)h-( a + /5 + l)n-j' 

Theorem 3.1. For m,n,s,t = 0,1, ... and (3 > —1, we have the orthogonality relation 

Z^ n {z,z)Z^ t \z,z)(zz) 0 e~ zz rdrde = tt — 1 ^ mAt- 

This follows from Theorem 2.1 in the case (3.2) but for completeness we give a direct proof. 
Alternate proof. When m > n we consider the integral 

I(m,n,a,b) := f Z^ n (z,z)(z a z b )r 20+1 e~ r2 drd6, 


for a < m,b < n. Therefore 

jY„, u\ — X (^1)™ f U \ (/^ + 1 )m ( ^ 

I(m,n,a,b) .- 27r<5 6+m , 0+n ^ Z. (jfeJ (£ + i) m _ fc ( 

n / \ / -i \n—k 


— l) fc / r . a + b + m +«-2fc+2/3+l e -r 2 ( ^ r 

0 


* (/3 + l)m \ 

— ^Ofo+m.a+n-j- 


n\ (-1) 


re! ^ VV (/3 + 1 )m-k 


r(l + (3 — k + (a + b + m + n)/2) 


= I rW,a + ..^tik £ (f) 7 tFZl f(l + p - k + a + n) 


n '- ^ W (P + f)m—fe 


= 7r(5b +m , n+n + ^ ^ + ^ m 2 Fi(-n, l + a + /3;/3 + l + m-n;l) 


— vrd^_|_ m]a _|_ n 


(/3 + l)m—n n! 

r(l + f3 + a) (f3 + l) m (m — n — a) n 
(f3 + l) m - n n\ (/3 + 1 + m - n) r 


Since m — n — a = —b we that the integral vanishes for b < n. If b = n then a = m and we 
conclude that 

f Zffl n (z, z)(z a z b )r 20+1 e~ r2 drdd = ttT(1 + f3 + m)5 a , m 5b, n , 

J R 2 

and the theorem now follows. ■ 

In view of (3.2) the recurrence relations (2.7), (2.5) and (2.12) become 

ziZ^ n {zi, z 2 ) = Z^l hn (zx,z 2 ) - Z^ n _ 1 (zi,z 2 ), 

*2 Z m+l,n(*l> Z 2) = ~( n + V Z m+l,n+l( z l> Z 2> + (P + m + 1 )Z^ n (z 1 ,Z 2 ), 

and 

(P + m + n + 1 - 22 ) 

= (« + i^m+hn+l^l’^) + (m + /3)Z^. lin _i(zi,«2) 
respectively, where m > n. We now discuss differential recurrence relations. 


(3.3) 

(3.4) 
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Theorem 3.2. For m > n, the polynomials {Zm,n(zi, £ 2 )} satisfy the differential recurrence 
relations 


S z 2 Z$ n (z 1 , z 2 ) = -z 2 Z^ ) n _ 1 (z 1 , z 2 ), (3.5) 

z 2 ) = nZ$ n (z!,z 2 ) - + m)Z < £ ) _ l n _ 1 (z 1 , z 2 ), (3.6) 

S z!Z$ n (zi, z 2 ) = (m- n)Z$ n (z 1 , z 2 ) - z 2 Z^ n _i(zi, z 2 ), (3.7) 

^i Z m)i( Zl ’ z2 ) = mZ mU Zl ’ Z2 ) - ( m + P) Z m-l,n-l( z hZ2), (3.8) 

($*1 - S Z2 )Z$ n (z 1 , z 2 ) = (m - n)Z^ n (z lf z 2 ), (3.9) 

(n5 Zl -m5 Z2 )Z$ n (z 1 ,z 2 ) = (m - n)(m + P)Z i £\ n _ 1 {z l , z 2 ), (3.10) 


where 8 Zl = z\d Zl , 5 Z2 = z 2 d Z2 . Moreover they have the operational representation 


Z m]n( z i’ 22 ) = ^exp(-9 Z i9 aa )zf +m ^ 

and the Rodrigues type representation 

45, = i 2 r' 5 e»»aj( i: r + ' , e—®) 

■( 2l j 2 rV>«s”((2 lZ2 )' , as( e -“‘«)). 


(3.11) 


n\ 

(-! r 


n! 


(3.12) 


Proof. Observe that from (3.1) we obtain 

(-I)”'" (/3 + l)m „m—k 


d z 2 Z$ n (zi,z 2 ) = 


k =0 


fc!(n - A:)! (/? + l) m _ fe 


^- fc (n - k)zkf 


n—k—1 


Therefore we have 


71—1 


d Z^ (zi Zo) = - — —_ (/3 + 1 )m rn-fc n-fc-1 _ _ 2 '(/ 3 ) z W 

22 m ’ n( i’ ^ ^ fc!(n-l-fc)!(/9 + l) m _ fc 1 ~ 2 i’ 2j 


/c=0 


and 


^*2-(^1.22) = Z~Y 


(-1)" (p + l) m m—k n—k 

Z-t Zrx 


z 2 ^ 0 k\(n-k)\(P + l)m-k 1 

(—l) n ~ fc (/?+l)m „m-fc n-fc-1 

" (k - l)!(n - A:)! {P + l) TO -jt 1 2 


m—l—j n—j —2 


- 22 m,n (h’ + Z. j!(n _ 1 _ i}! (/3 + 

n ^ , (P + m) ^ (-l)"- 1 -^ (0 + l) m - l2 r 

= ,2) - Y j !( n — 1 — j)!- (/3+ 


£2 “) j!( n — 1 — j)! (/3 + l)m-l-j 


_ n -7(A) r \ (/3 + m) (P) 

^m,n\Z 15 -^ 2 ) ^m —1 ,n—l (^ 1 ? ^ 2 ) 

^2 ^2 


The relationships (3.7) and (3.8) are proved similarly. Formulas (3.11) and (3.12) follow by 
direct calculation using (3.1). The first-order partial differential equation (3.9) can be obtained 
from either (3.5) and (3.7) or (3.6) and (3.8), (3.10) is proved similarly. ■ 
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One would have expected the Rodrigues formula for Z^ n (z i, z 2 ) to be a constant multiple of 

which is obviously false except in the 219-Hermite case P = 0. 

Note that (3.4) follows from (3.12) by writing d r f as dff 1 d zi . 

We shall use the notation 

wp(ziz 2 ) = {z x z 2 )^ exp(-z x z 2 ). 

It is clear that (3.12) implies the differentiation formulas 

e zlZ 2 d Z 2 (e ~ zlZ2 Z^ n ( Zl , z 2 )) = -ziZ^+ l \zi, z 2 ) = ~Z^l hn ( Zl , z 2 ) 

and 

wp(ziz 2 ) dzi ( W P( ZlZ 2 } Z ™U Zl ’ Z ^) = ( n + l ) Z m,n+l( z l>z 2 ) (3.13) 

for m > n. It also clear that 


d Zj wp{ziz 2 ) = (P~ ziz 2 )wp(ziZ2), j = 1,2 

and (3.3) which imply the relationships 

6z 2 {™p(z 1 z 2 )Z$ n (z 1 ,z 2 )) _ (/3) 

l \ P Z m,n\ z li Z 2 ) n (~l j ~2); 

wp{ziz 2 ) 


(3.14) 


and 


S Zl {wp(z 1 z 2 )Z^ n (z 1 ,z 2 )) _ (0) , 7 C8) , 

yjp(^Z\Z2) (p + 7T1 Tl)Z/ m n yZi^ Z 2 ) ,nV^l’ 


This can be seen by first applying (3.14), then applying (3.3), for example, 
1 


wp(z x z 2 ) 

P - ZlZ 2 


d Z2 {wp{z l z 2 )Z^ n {zi,z 2 )) 


Z $n( z hZ 2 ) - Z^ n _ 1 (z 1 ,Z 2 ) = ^ r Z$ n {z l ,Z 2 ) - Z^l l n ( Zl ,Z 2 ). 


p 


d/5) 


^2 ’ 22 

Theorem 3.3. The polynomials {Zm}n(zi, z 2 )} satisfy the second-order partial differential equa¬ 
tion 

'P - ZlZ 2 


ddz i &d Z 2 f + 


21 


ddz 2 .f = -nf, 


(3.15) 


for all m > n. 

Proof. From equations (3.13), (3.14) and (3.5) we get 
d Zl wp(ziz 2 )Z$ n (zi, z 2 ) + wp(z 1 z 2 )d Zl Z^ n (z 1 , z 2 ) 


wp{z x z 2 ) 

= Z2 ) Z m.n(2l, 22) + d Zl Z$ n (z X , Z 2 ) = (n + l)Z^ n+1 (zi, Z 2 ) 


and 

dz 2 { (£ - 22 ^) PS,( 21 ,22 ) | + d Zl d Z2 Z$ n (zi, Z 2 ) 

= (n + l) 9 * a ^ i ) n+ i( 2 1 ,22) = -(n + l)Z^ n (2 1 , z 2 ), 
which simplifies to the desired second-order partial differential equation. ■ 
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It is clear that the differential property (3.15) can be written in the form 


d_ 

dz\ 


d 


wp(z 1 z 2 )—Z^ n (z 1 ,z 2 ) = -nZ$ n (zi, z 2 ) 


(3.16) 


Note that (3.15) and (3.9) indicate that the polynomials {Zm} n (zi, z 2 )} are simultaneous 
eigenfunctions of the operators on their respective left-hand sides. The differential operators 
d dzi dg Z2 + ( /3 ~" 1 122 ^ dg z2 and z\dg zx — z 2 dg Z2 indeed commute as can be directly verified. 

Let w(x,y) = (zz)^e~ zz = r 2 ^ +1 e~ r ~drdO. It is clear that the bilinear functional 


(/. a) 


f (x, y)g(x, y)w(x , y)dxdy 


defines a semi-inner product on the linear space of w(x,y)dxdy measurable functions f(x,y ) 
such that f R2 \f(x,y)\ 2 w(x, y)dxdy < oo. Since a polynomial p(x,y) is also a polynomial in r 
with trigonometric polynomials as coefficients, hence it is a continuous function in ( r,6 ). Then 
f R 2 \p{x, y)\ 2 w(x, y)dxdy = 0 implies that p(x,y) = 0, consequently the bilinear functional is 
positive definite on the subspace of all the polynomials p(x, y). In the subsequent discussion we 
only need this inner product space of all the polynomials. 

Theorem 3.4. We have the adjoint relations 

(d z y = -dz - ^ + Z, {d z y = -d z -^- + z. 

The proof is straightforward calculus exercise. 

This allows us to write (3.15) or (3.16) in the selfadjoint form 

m)*d z )f = nf. 


Moreover we also have the following result. 

Theorem 3.5. The operator A := (■ d z )*d z is positive in the sense that (Af , f) > 0 with = if 
and only if d z f = 0, that is f depends only on z. 

This indicates that 

1 d Z m]n ( Z » z) = (n+ l)Zff n (z, z). 


w(x, y) dz 


This is the adjoint of (3.5). 

Ismail and Zeng [20] established the connection relations stated in the next theorem. 
Theorem 3.6. We have the connection relation 


Z mjn( z i. z ?) = 


(P - l)i 


(- 1 YZ™_j jn _j{ Zl , z 2 ), m > n. 


i =o 


For m > n, we have the special cases 


n!z m,n(*i» z 2 ) = 'jr f™) (/3) j (-iyn m - j ^ j (z 1 ,z 2 ), 

d—n / 


j =0 

n 


H m ,n(zi,z 2 ) = n\ ^2 Z^l jtn _yzi, z 2 ). 


(h) 


3=0 
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The Laguerre differential equation is [28, 31] 
xy" + (1 + a — x)y + ny = 0 . 


Therefore Theorem 2.10 shows that 

0 ) 

^- + (1 + /3 + m — n — Z 1 Z 2 ) 


Z2~ 


a2 y{0) 
u ^m.n 


dz\ 


v £Jm,n 

dzo 


+ nziZ^ n = 0 . 


Ismail and Zeng [20] gave the generating function 


£ 


(m — n) 


m>n> 0 
this can be seen from 


] z m, kOi,z 2 ) = (1 + uv) p 1 exp , 


g Q (x,z) = f2L<?+V(x)z n = ( 7^ g )a+/+l ' 


and ( 2 . 21 ), 

E 


m>n>0 


(m — n)\ 


Z mU z ^ z 2) = 


r uz 1 -z 1 z 2 uv \ 

fcX P l 1 — uv I 

(1 - uv) p+1 


pxn / _ Z1Z2UV \ 

E u m v“ZWfr,Z2)= „ l, 1 -™' - T. 

^ ’ {l-uv)P(l-Ziu-uv) 


m>n> 0 


(3.17) 


It is clear that the generating function (3.17) implies the identity 


Z m,t T+1) (2l +23,22 + 24) 


£ 

m>j>k> 0 


(^2)EEn-fc(~3, Z4) 
(j — k)\(m — n — j + A;)! 


This is an analogue of the convolution identity 


47' 5+1 >(+ + !/) =E 0“\x)Lf>(ji). 

k =0 

Al-Salam and Chihara characterized all 1 -D orthogonal polynomials satisfying convolution for¬ 
mulas in [1], They discovered the Al-Salam-Chihara polynomials through this characteriza¬ 
tion [17]. It will be interesting to solve the corresponding 2 -D characterization problem. 


4 The polynomials z 2 )} 

For a > 0, /?, 7 > — 1 Ismail and Zeng [20] introduced the polynomials </> n (r, a) 
that is 


p(« + 7A)(i_ 2r ), 


r, a = 


(a + 7 + 1 ) n ^2 


n—k 


k =0 


(a + /3 + 7 + n + 1 ) n _ fc (-r) 
k\(n - k)\{a + 7 + l) n -k 


(4.1) 


They satisfy the following orthogonality relation 
cl 


10 


<t>m(r, a) 4 >n(r, a)u a+1 (1 - u) p du = ( n (a + 7 , P) 5 m ,r 
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where 


r(q + 7 + n + l)r (/3 + n + 1) 

n a n'.T(a + (3 + 7 + n + l)(a + j3 + 7 + 2n + 1)' 

We define the two variable polynomial z 2 ) by 

7 ,^, 7 )/ N (m +/3+ 7 + l) n -fc( 7 +l) m z^ k (-z 2 ) n ~ k 

M m , n (*1, ,2) = X. -- 


for m> n and 


} (*i> - 2 ) = (* 2 , *1), m < n. 

Then 

^l A ^n +1) (-1^2) = ^m+l,n( z l^ 2 ). 

Clearly, Mm,n\zi,z 2 ) satisfy the orthogonality relation 



_ T (7 + m + l)T(/3 + n + 1 ) 

n\T((3 + 7 + m + l)(/3 + 7 + m + n + 1 ) m,p n,q ’ 

for m > n as it is stated in [20]. The disk polynomials defined in [7] can be expressed as 

p lL,n( z ) = *)’ m > n. 

From (4.1) we get 


(m + /3 + 7 +l) n _ fc ( 7 +l) m (~l) n 

Ct( "'“ ) = -t!(»-m( 7 + l) m . t -’ 

where a = m — n with m> n. Then 

(/3 + 7 + m + n + 2)z 2 M^^ n (z 1 , z 2 ) 

= (7 + m + l)MVtf(z u Z 2 ) - (n + z 2 ), 

(/? + 7 + m + n + z 2 ) 

= (0 + 7 + m + l)M, ( 7 f + 7 1 ) .„(zi, z 2 ) - (/? + z 2 ) 


and 


(c n - ziz 2 )M^\z t ,z 2 ) = a n M^l n+1 ( Zl , z 2 ) + bnM^^Zt, z 2 ), 


(4.2) 


(n + l)(/3 + 7 + m + 1 ) 

(/? + 7 + m + n + l)(/3 + 7 + m + n + 2 )’ 
(7 + m)(/3 + n) 

(/3 + 7 + m + n)(/3 + 7 + m + n + l)’ 


where 
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(n + l)(7 + m + l) 71(7 + m) 
n [3 + 7 + m + n + 2 j 3 + 7 + m + n 

It must be noted that (4.2) is essentially the three term recurrence relation for Jacobi polyno¬ 
mials [31]. 

It is clear that 

A n 5 l<t> n (x, a) + B n 8 x <j) n (x , a) + C n (j> n (x , a) = 0, 

where 

A n = l-x, B n = a + 7 — x(a + /3 + 7 + 1), C n = m(a + /3 + 7 + n + 1). 

Then for m > n, Mm’n\zi, 2 2 ) satisfy the following second-order partial differential equation 

(1 - zi22)6 z 2 /(2i, 22) + {m - n + 7 - 2122O8 + 7 + m - n + l)}S Z 2 f(zi, z 2 ) 

+ ziz 2 n (/3 + 7 + m + l)f(zi, 22) = 0 . 

An equivalent form is 

n \ d 2 f(z!,z 2 ) , ri , , /0 , , « , ^df(zi,z 2 ) 

(1 - ziz 2 )z 2 -—^-hl + m- n + 7- z 1 z 2 (2 + 7 + p + m- n)\ --- 

OZ2 oz 2 

+ zin(m + /3 + 7 + l)/(,zi, z 2 ) = 0 . 

Similarly we establish 

n ^ 9 2 f{z 1 ,Z 1 ) , m<9/(2i,Z 2 ) 

(1 - 21^2)21-—j- Hl-m + n + 7- ziz 2 (2 + 7 + p - m + n)\ --- 

OZ2 oz 2 

+ z 2 m(n + /3 + 7 + l)/(zi, 22) = 0 . 

Theorem 4 . 1 . For m > n > 0 , /3 > — 1 , 7 > — 1 we have 

■^ M m,n ] ( z 1. 22) = -(^ + 7 + m + 1)M^+ 1 ^( 2 i,^ 2 ), 

$z 2 M^fi\z 1 ,z 2 ) = nM^\zi,z 2 ) - (7 + m)M^ f 1 1 ^l 1 (2i, z 2 ) 

= (/3 + 7 + m + 1 )(M^+ 1 > 7 ) (21,22) - M^\ Zl ,z 2 )), 

M m,n\ z ^ z 2) = mMV$\zi,z 2 ) - (7 + « 2 ) 

= (/3 + 7 + m + 1)M^+ 1)7 )(2:1,22) - (/3 + 7 + n + 1 )M^ (zi, z 2 ), 

(P + ^ + m + l)S Zl M^\zi,z 2 ) - (P + 'y + n + l) 6 Za M}Pg\z 1 ,z 2 ) 

= {m - n)(0 + 7 + m + 1 )M^+ 1 ’^(z 1 , z 2 ), 

(S Z1 - S Z 2 )M^\zi,z 2 ) = (m-n)M^\zi,z 2 ), 

(/3 + 'y + rn + l) 6 Zl M$g\zi,z 2 ) - (/3 + 7 + n + 1 ) 5 Z 2 M^ (z 1: z 2 ) 

= (m- n )(/3 + 7 + m + l)M^ 1 , ' y \zi, z 2 ). 

Proof. The first equation easily follows from the definition. From the definition iorMm,n\zi,z 2 ) 
we get 

r t\,t(0,i)( _ \ (/3 + 7 + 771 + l)n-fc(7 + l)m2] 71 k {—Z 2 ) n k (n — k) 

6 „M m , n ( ZI , 22 ) = ^-fc!( n -*)!( 7+ l) m - t - 
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(/3 + 7 + 171 + l)ra-fc+l(7 + k {—Z2) n k 

“ A:!(n- fc)!(7 + l) m _fc 


-(^ + 7 + m + l)J] 

fc =0 


(/3 + 7 + m + l)n-fc(7 + l)m*™ k {—Z2) n k 
k\(n - k)l (7 + l) m _fc 


= (/3 + 7 + m + l)M^+ 1 ’' y \z 1 ,z 2 ) - (/3 + 7 + m + l)M^\z u z 2 ) 


x T\A0n)( (P + 7 + m + l)n-fc(7 + k (—z 2 ) n k (n — k) 

S„M m-n (*„ z 2 ) = Y-*!(„ - fc)!(T, + l )m _ t - 

_ (/3 + 7 + m + l) n -fc (7 + l) m z™ k (—z 2 ) n k 

H “ fc!(n - fe )!(7 + l)m-fe 

_ y- (/3 + 7 + m + l)„._fc(7 + k {—z 2 ) n k 

“ (fc - l)!(n - fc)!(7 + l) m _ fe 

= n- M m,n\zi,z 2 ) - (7 + z 2 ). 

Similarly, 


y 7 \f(/ 3 , 7 V~ „ \ 'sy (/5 + 7 + m + l)n-fc(7 + k (—z 2 ) n k (m-k ) 

(ZI ’ Z2) = E-*:!(»-fc)l(7 + D™-*- 

(/3 + 7 + m + l) n _/i;(7 + l) m z] Tl fc (—2:2) n k 


= m 


E 

fc =0 
n—1 

E 


fc!( n “ k)!(j + l) m _ fc 


\n—1—/c 


fc =0 


(/3 + 7 + m + l)n—1—fc(7 + l)m~i n ( — Z 2 ) r 

k\(n - 1 - fe)!(7 + l) m -i_fc 

= ™MW$(zi,Z 2 ) - (7 + *2)> 

x \ (P + 7 + m + l)n-fc+i(7 + l)mZ™ k {—z 2 ) n ~ k 

4 1 <^(, 1 ,z 2 ) = g- fc! ( w - fe )!(7 + i) m _ fc -(/3 + 7 + n + l) 

_ (/^ + 7 + m + l) n -fc(7 + l) m 0™ k (—z 2 ) n k 


k =0 


fc!(n — fe )!(7 + l) m -fc 
= (/3 + 7 + m + 1)M 1 ^+ 1 ’ 7 ) («i,« 2 ) - (/3 + 7 + ra + 1 )M^ ) (zi,« 2 ). 
Ismail and Zeng [20] used the generating function 
x 9 a+/3 

n=0 


R(l + ^ + R)/ 3 (l-^ + R) c 


(4.3) 


where 

R = \J (1 — z) 2 + 4 x 2 , 

[28, Section 140] to establish the following generating functions: 

00 9/3+7 

y: Z2)u m V n = - (1 + uv + p) -/3 (1 -uv + p) -7 

m,n=0 ^ 
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+ 


and 




U m V n 2^ +7 


m — n)! p 


1 — 2^irt/(l + p — uv ) 1 — 2 vz 2/(1 + p — uv) 

(1 + uv + p)~ 13 (1 — uv + p)~ 7 exp 


- 1 


2zi u 


m>n> 0 

where 

p = (1 — 2 ieu(l — 2Z1Z2) + u 2 v 2 ) 1 ^ 2 
They also verified the limiting relation 

lim M^\ Zl ,z 2 /P) = z?- n LX +m - n \ Zl Z 2 ) = (-1 rZ$ n (z 1 ,z 2 ). 

p—too 

One can use the generating function (4.3) to establish the generating relation 

u m v n M^\ Zll Z 2 ) = 2/3+7(1 + UV + P)_/3(1 " UV + p)l ~ 7 

m>ri> 0 

with p as in (4.4). 


1 — uv + p 


(4.4) 


p(l — uv — 2uz\ + p) 


5 The polynomials 1<?)} 

For a > — 1, the moment problem associated with g-Laguerre polynomials 

” 1 ,q) [q (q\ q)k(q,q a+1 \q)r 


is indeterminate. It is well-known that 


)n—k 


TOO 

/ L m\x; q)L { n\x-, q)x a dp(x\q ) = C n{a)5„ 
Jo 

where p(x\q), Cn(«) are given by 

Cy. /v* 

dp(x\q)= . . 


(~x;q )oo ’ 


(g a ;g)oo (g a+1 ;g) ? 

(g;g)oo q n (g-,q)n 
xd(x — cq k ) 

k =-oo ’ 

(g,- C(? Q+1 ,-^;g) oo (g Q+1 ;g) ? 


Cn(a) = r(—a)r(a + 1) 

OO 

dp{x\q) = 

Cn{oc) = 


(g a + 1 ,-c,-2;g) c 


—a:—1 


(g; g)ng^ 


and 


dp(x|g) = 

Cn(®) — d mn I 
Jo 


(-**>-&; 9)oc ds 


(-x,-Xq c x,- J 3 E --q)J 

x a ~ c {-\ s,~^jg^rfs 
(-x,-A q c x,-^q)J 


(5.1) 


(5.2) 


c, A > 0, 


(5.3) 
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where c, A > 0. Following the procedure outline in Section 2 we let 


42i(-i,-2k) = 


z™ n Ln +m n \ziz 2 -,q), m > n, 


z$h( z 2,zi\q). 


m < n, 


then we have 


n _ V q ^ +m ~ k){n ~ k) (Q ni ^)mzT~ k z. 

^ h (- 1 ) n_fc (9; ?)fc( 9 ; q)n-k(qP +1 ;q)m-k 


—k„n—k 
2 




and 


W? Q.)n 

for m > n. Applying Theorem 2.1 we obtain the following result 

Theorem 5.1. For m,n, s,t = 0,1,... and f3 > — 1 we have the following orthogonality relations 

[ Zm] n {ZiZ\q)zft( z ,z\q)r 2 Pdii(r 2 \q)dQ = ( mA n(.\m - n\ + /3)5 m , s d n ,t, 

Jr 2 

where d/a(x\q), C,n{ot) may be any pairs (5.1), (5.2) or (5.3). 

Applying formulas (2.5) and (2.7) in the case 

tqOL+l. g\ q(a+n—k)(n—k)r^n—k 

c k (n,a) = - ’ , -—TA- ( 5 - 4 ) 

{q-,q)k{q,q a+1 -,q) n _ k 

we obtain the recurrences stated as a theorem below. 

Theorem 5.2. For (3 > — 1 and m>n we have 

q {m+1+l3) Z 2 ) = -(1 - Q n+1 )z^l 1>n+1 (zi,z 2 ) + (1 - q m+ P +1 )zffl n (z 1 ,z 2 ), 

q n ZlZ^ n+1 (zi,Z 2 ) = z£l lin (zi, Z 2 ) - i ffl n _l( 2 l,Z2) 

and 


(-i y 


(1 + q{l -q n - q m+ P) - z lZ2 q^ m+n+1 )z^ n ( Zl , z 2 ) 

= (1 - q 1+n )z^_ l n+1 (zi, z 2 ) + q(l - g /3+m )^! lin _ 1 (AA 2 ). 
The q-Laguerre polynomials satisfy the following ^-difference operators 

D q>x L^\x-q) = -L [ ^\qz;q), 

i _ qTi+I 

D q ^{w a (x-q)L ( f x \x;q)} = w a _i(x; q)L^\z- q), 

w a {x-q)L^\x-,q) = ^ ^ D™ x {w a+n (x; q)}, 

\Qi Q)n 


and 


wg(x-q) 


x^ 

(~x\q) oo’ 


(5.5) 

(5.6) 
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Theorem 5 . 3 . For /3 > — 1 we have 

D q ,zi{zi~ m z$ n (z 1 ,z 2 )) _ q 13 z 2 (,3) 


q- i-z i 


z^ n _ 1 (qz 1 ,z 2 ), 

1 Zl 


^q,z 2 ( z m.,n (■‘'1) ^ 2 )) _ ^ ^m,n— 1 > ^ 2 ) 1 


9 


z m}n( z ^-i ^ 2 ) 


(1 - ^) m+n 0- 2i {z|0™ 2 {^ /3 (z 1 z 2 )}} 


(—l) m ((/; q)n'wp(ziz 2 ) 


Z2D q!Z2 {wp(ziZ 2 )z$ n (Zi,Z2)} _ 1 ~qP 


,tw /v. _ ffjfLJP) ( z z ) 

z m,n\ z 1? _ 2 m+l,nv 2 l’ ^2,/ 


wp{ziz 2 ) l- q 

for m > n and 

Dq iZl {wp(z 1 z 2 -,q)z$ n (zi,z 2 )} _ 1 - q n+1 


"m,n+l(^l’ ^ 2 ); 


W/3(ziz 2 ;q) l-q 

D q,z 2 {z?~ n Wi3(ziZ 2 -,q)z$n{zi,z 2 )} _ 1 - q n+1 Zi 

z™~ n wp(zxz 2 -,q) l-q z 2 

for m> n. 

Proof. First we observe that 


y (/J) 

'"m,n+l 


( 21 , 22 ) 


(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 


1 


D q , x {x Pwp(x)} = - q —^ X Pw p( X ) 


to get 


z 2 +rn D^z 2 { z 2^ w p(ziZ 2 )} = 


Wp+m{ziZ 2 ) 

(q- l) m 


and 


or 


P+m n m { p __ u -| _ D g, Zl { W P+m(ziZ 2 )} z r i 


D q,zA Z 2 D ™z 2 { Z 2 Wp( z l z 2 )}} = 


(q-iy 


( 9 - 1 ) 


m D q,ziz 2 { w P+m( z l2 2 )}, 


^ 2l {4 +m " n ^ 2 {^^(^1^2)}} 

_ _ (~ 1 )"‘(g; g)n n„(/3) / \ 

— ( g _ ])m — (1 _ g)m+ra ~ 2 z m,nOl) Z 2j, 

which is (5.9), (5.7) and (5.8) are obtained from (5.5), and (5.11) and (5.12) are derived 
from (5.6). 

From (5.9) we get 

(1 _ \m+n 

D q,z 2 { Wp( z l z 2 ) (^1, *2) } = (_ 1 )m( g . g ) ^9,*! { ^9 ,*2 { ^2 ^za { Z 2^ W P i Z l' Z 2 ) ) ) ) 

(l_ g) m+n-l(i_^) nn f j3 _ lnm 

q) n ^q,zi\ Z 2 ^q,z 2 \ Z 2 w p(~'l z 2)}} 


+ 


(1 - q) m + V nn 

(-i) ro (g;g). 


K^t 1 K%(^ 2 )}}, 


which gives (5.10). 
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The g-Laguerre polynomials satisfy the following second-order difference equation 
-<?"(! - q n )xy{x ) = q a (l + x)y(qx) - (l + q a (l + x))y(x) + y(q~ 1 x), 


or 


(1 + qx)0 2 q . x y(x) + 


1 — q a — 2g“ +i .T + q a + n + l x n / N | qx( 1 — g") 


g Q (l-g) 


-O q , x y{x) + 


(1-9) 2 


y(x) = 0. 


where y(x) = L^\x;q), 6 QtZ = zD q>z . 


Theorem 5.4. For m > n the function f = ZmL( z \ , z 2 ) satisfies the second-order q-difference 
equations 


{l + qziz 2 ) 9 \, 2 f + 


- 1 - (2 - q n )qz 1 z 2o f , g*iz 2 (l - q n ) f _ n 

~dq,z 2 J + n x 2 J ~ U 


1 ~q 


(1 ~qf 


and 


{l+qziz 2 )e q ^ z J - 


2 i _ qP+m-n + ( 2 _ qh+m+l) qZl z 2 q(l-qh+ m ) Zl z 2 


1 - q 


1,21 f + 


(i - qf 


f = o. 


Starting with the little g-Laguerre [23] or Wall polynomials [5] 

pn(x . q « lq) = y iMU -*- t— k 

“ ( q-,q)k(q,q a+1 ',q)r 


)n—k 


we define the 2D polynomials Wm] n {z\, z 2 \q) through 


Z 2 \q) (yg) 


* 5 " n Pn(-2l22;^ +m n \q), m> 


n. 


K Wn,rn(z 2 ,zi\q), m < n, 

for f3 > — 1 and m,n = 0,1,.... Then for m > n find that 

i, z 2\P) = w£l ln { Zl ,Z2\P), 


w, 


( 0 ) 


\{zi,z 2 \q) = ^ 


k =0 


(g;g)n(g^ +1 ;g)m-n-Sr~ fc ^~ fc g~ ( * ) ~ (S) 

(-i) n_fc (g; g)fc(g; g)n-fc(g /?+1 ; 


The orthogonality relation for the little g-Laguerre polynomials is 

OO 

Yl q (a+1)k (9 k+1 -,q)^p m {q k ;q a \q)p n {q k -,q a \q) 


k =0 


( q;q)ocq (a+1)n (q;q)nSn 


(q a+1 ;q)oc(q a+1 ;q) n 

This leads to the orthogonality relation {wm,n(zi, z 2 \q)}™ n=0 in the form 


m,n = 0,1,.... 

. 08 ) 


Theorem 5.5. For m,n, s,t = 0,1,... and (3 > —1 we have the following orthogonality relations 


w m]n( z i Aq) w s,t( z ’ z\q)r 2l3 dy(r 2 \q)d0 = CmAn{\m - n\ + /3)5 m , s 5 n , t , 
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where z = re ld , 

OO 

dn( r \Q) = ^ ~2x(xq-,q) 00 8(x - q k ), 

and 


k =o 


C nipt) — 


( q-,q)ocq ia+l)n (q;q) J 


(■ q a+1 ',q)oo(q a+1 ;q)n 
The function y(x) = p n {x', q a \q) satisfies the (/-difference equation 

q a y(q 2 z) - (l + q a - V~ n )y(<?T> + (1 - qz)y(z) = 0 
or equivalently 

q nl + q° 


q a+n - l 0i z y{z) + 


-,n- 1 _j_ „a+n-l _ ^ z(l - q n ) 


— Z 


8q,zy(z) 


' ' 1 — g ' ' (1 — g) 2 


y{z) = o, 


where 9 q , z = zD q ^ z . This leads to the following theorem. 

Theorem 5.6. For (3 > —1, m>n the polynomial Wm) n (z\, z 2 \q), satisfies the q-partial diffe¬ 
rence equations 

0+m—lo2 , q n ~ l - q 0+m ^ - z lZ2 (l-q n ) Zl Z2 

Q Vq,z 2 J H -i-7- J H-71- 


i - q 


(i - qf 


and 


„ no 2 , q n - q l3+m + qziz 2 Q f , g(l ~ q^ +m )z lZ2 f _ n 

Q “q,ziJ 1 „ “q,ZlJ + 9 J — U. 


1 - q 


From 


D q ^xPn{^iq q ) — 


g 1_n (l - g n ) 


(1 - qf 


-p n ^l(x;q a+1 \q), 


(l-q)(l-q a + 1 ) 1 
D q -\ x {w(x-q a )p n (x-,q a \q)} _ {l - q a )p n+l (x] q a ~ l \q) 


w{x ; a — 1) 


<? a_1 (i - q) 


(5.13) 


and 


q na H 2)(1 - g) n 

w(a;;a|g> n (x;g a |g) =-- D”-^ x w(x; a + n\q), 


(q a+1 -,q)r 


where 


w(x-, a ) = (qx; q)ooX a 
we get the following relations: 

Theorem 5.7. For ft > — 1 and m > n we have 

D q,zi{zi~ m 'W^ n {z 1 ,Z 2 \q)} _ z 2 ^“"(l “ l(«l> ^l?) 

z n-m - ( 1 -g)(l_ g /3+m-n+l) ’ 

9 1-n (l - 


£) 9,*a{ w m l n( 2; l> 2! 2|9)} = “ 


(1 - q)( 1 - qr/3+m-n+l) 
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and 


w(z 1 z 2 \P\q)w ( £] n (z 1 ,z 2 \q) = 


^m(n-l)+n(j3-l}-(j)^ _ ^m+n 
(-1 )m(^+m-n+ 1 ;? ) n 


For m > n we have 

D q-\z 1 { W ( Zl /^|g)Wm,l(^l,^ 2 |g)} _ (1 - g j3+m " n )'^!n, ) n-n(^l ; ^[g) 

w(2;iZ2;/5|g) g/9+m-n- 1 (l _ g ) 

^-h, 2 {^ ra_n ^(^i^;/ 3 k)^£U^i^ 2 k)} _ zi (1 - g / 3 +m - n )T0^ ) n+ 1 (2 1 ,z 2 |g) 


4 1_n w(2;i22;/3|g) 

^2Fg-i jZ2 {u;(2;i2;2;^|g)^m,U^i,^2k)} 
-221<?) 


^2 


q P+m-n -!(]_ _ g) 


1 - q 

Proof. Observe that 


1 ~q p (0), ^ ^ +1 " n 1 - q 8+m+1 09) , 

— ^/n(R^2k) - x _ x _ /3+m-n+l g 2W m +i |tt (^, ^Ig) 


n™ — v n D n 
Uq-l, Zl ~ Z 2 U a ~ 


q~ ± ,ziz 2 


and 


Z?W( Z 1 Z 2 ;/%) = (l - g- 1 ) m (z 1 z 2 )^ m _ 1 


then 


9 n(^+m)-("+ 1 )( 1 _ g ) 


W(^l22;/3|g)w^ 1 ) n («l>^2|9) = - ( g fl +m - n+1 . - D q-\ zl { z Tw{ Zl z 2 -, P\q)) 

q m {n-l)+n(p-l)-(™) m _ Jm+n 

- (_l)m( g /3+m-n+l ;g ) n 

X { (212 2 )' /3 W(Z12 2 ; /?|«)}}, 

? m(n-l)+n(/3-l)-(”)^ _ ^m+n 

(-l) m (g ^+ m - n+1 -,q) n 

* D n rlzi {(z 1 z 2 fD™ t ^{(ziz 2 )-l , w(z 1 zr,P\q)}}, 


w(zi Z2;0\g)w l ;P n {z i ,Z2\q) = 


which gives (5.14), (5.15) and (5.16) are obtained from (5.13) by direct computation. 
From (5.14) we have 


D q- 1 ,z 2 { W { Z ± Z 2\P\q)w { JSn(zi.,Z2\q)} = 


g m(ra-l)+ra(0-l)-(”)^ _ ^m+n 


(5.14) 

(5.15) 

(5.16) 


(5.17) 


(-1 )m(^+m-n+ 1 ;g ) n 

m(n-l)+n(/3-l)-(”)/-. _ \m+n 

{ -rr( q 0*r,J + uJ n 

g m(n-l)+re(/3-l)+l-0-(”)^ _ ^m+ra-l^ _ g /3) 

+ (-l)m(qh+m-n+l. q ^ n 

X ^2 _l£, 5-l, 2 l {(R^) / 3 F”i 1 ) 22 {(2l2 2 )' / 3 W;(2l22;/3k)}} 
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qP+l-n l _ q 0+m +1 

~ 1 _ q l _ qP+m -n + M Z ^ P\9) w m+l,n( z l’ 
^ Q ^ w(ziZ 2 -,P\q)w^Uzi,z 2 \q), 


z 2 (l - q) 
which gives (5.17). 

Applying (2.8)-(2.11) and (2.13) to 


c k (n,a) = 


(k-n)(n-\-k — 1) , 

(g;g) nq * (-1 ) n ~ k 

(q;q)k(q,q a+1 ;q)n-k 

we get the following recurrences: 

Theorem 5.8. For j3 > —1, m > n we have 

q) (0) , (0) , 

g n (1 _ q m-n+0+ 1) = W m!n . -2 \q) ~ <+l,n+l(*l> ^l?). 

= (1 - 9 m+ ^>Si, n (^i, 22 k) - g m -^ +1 (i - (f )<Li(^2k), 

and 

{q n + q m+P { 1 - 9 n - O - z 2 |?) 

= 9 n (l - 9 m+/3+1 )^i 1)n+1 (^i^2k) +g m+ ^(l -q n )w^\ n _M^2\q). 

6 The polynomials z 2 \q)} 

For a ,7 > —1, the Little g-Jacobi polynomials 

D ( x . a « a 7\ a ) = (gjg)» V (g a+wl ;g)n- fc g(^(-x)"- fc 
gGO “ (g;g)fe(g,9 a+1 ;g)n-fe 

satisfy the orthogonality relation [23] 

00 (~k{ot.-\- 1) 

J2 ( g 7+fc+";g) c -^’ 9 7 |?)Pn(? fc | 9°, <7 7 |</) 


k=0 


{q, q a+ ' y+n+ 1 ;q)oo(q', 9)n</ l(a+1) <W 


{q a+1 , q ry+n+1 ;q)oo{q a+1 ] q)n( 1 - g«+ 7 + 2 n+l)' 
This leads to the choice 




/c=0 


Cn(o,7k) = 


(q, q a+ ' y+n+ 1 ;q)oo(q', q) n q n( - a+1) 


(q a + 1 ,q'y+ n + 1 -q) 00 (q°‘+ 1 -q)n(l - q^+i+^n+iy 
Following our general construction we define the polynomials 

,, 09 , 7 )/ | \ j z^~ n p n (z 1 z 2 ,qf s+m - n ,q^\q), m> n, 


( 6 . 1 ) 
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Then we have the explicit form 

M {Pn)( | \ _ (g;g)n(g^ +1 ;g)m-n (g /j+7+1 ; q)m+n-kq^ Z^ k (-Z 2 ) n ~ fc 

m ’ n 1 ’ 19 (q^-,q)mq® h & ^ 9)n-k(q^ +1 \ q)m-k ’ 

z l M m,n 1 ’ l) (zi,z 2 \q) = M%£l n (zi, z 2 \q) 
for m > n. Applying Theorem 2.1 we obtain: 

Theorem 6.1. For m,n,s,t = 0,1,... and (3 ,7 > —1 we have the following orthogonality 
relation 


/ M^Q\z,z\q)M^ t ,1 \z,z\q)r 2 ^dfj,(r 2 ;'y\q)d0 = CmAn(\m - n\ +/3,j\q)S m , s S ntt , 
Jr 2 

where dfj,(x\ j\q), £ n (o:, 7 | q) are given in ( 6 . 1 ). 

Next we apply formulas (2.8)—(2.11) and (2.13) to 

{q-,q)n{q a ^ +n+1 -,q)n-k{-l) n - k 


Ck{n,a ) = 


(n\ /k\ 1 


(q]q)k{q-,q a+1 -,q)n-kq^ 2 > ( ' 2j 
and obtain the following recurrences: 

Theorem 6.2. For /3, 7 > —1 and m > n we have 

)!■(! - q fi+m-n+V) ZlM ^+l\n( Z1 ’ ^ = -M m+h+l( 2 l> «l«) + 

= (1 I 

f-wd-fid-n ,^. 7 ),, .. 

(7 _ g/3+m-n+lVj _ q/3+^+m+n+l j m,n— lv^l) 

The polynomials {p n (x; 9 a , Q^lq)} satisfy the following second-order difference equation 
(l + q a - ( q l ~ n + q n+a+1+2 )x)y(xq) = q a { 1 - q 1+2 x)y{q 2 x) + (1 - qx)y(x) 


or 


(1 - q 1+2 x)0l x y{x) + 


q a — 1 — q 1 71 a x — q’ y+2 (2 — q n ^x 


1 -q 

1 + q 1 ~ a x — (g 1_n_a + q'y +n+2 )x 


0 q ,x y{x) 


y{x) = 0 , 


(1 - v ) 2 

where y(x) = p n (x ; q a , <? 7 |g), 9 q , z = zD q , z . 

Theorem 6.3. For m> n and /?, 7 > — l, the function f = z 2 ) satisfies 


q n—m—[3 _ 1 _ q l-m-p ZlZ2 _ g7 +2(2 _ c f) ZlZ2 ^ 


° q ’ Z2f+ l (l-g)(l- 97 + 2 ^) 

1 + (q n + l ~ m -P - ql-m-P _ ql +n+2^ ZiZ2 


9q Z 2 f 


(1 - q) 2 { 1 - q^ +2 ziz 2 ) 


f = 0 


and 


3 2 ( 21^2 (2g 7+2 - q 1 " n - g™+d+7+2) + q m-n+P _ x 


<W + 


(1 - g)(l - 9 7 + 2 2 iZ 2 ) 


^ 1 / 


(q 1 ~ n _ gl+m-n+/3 _j_ ^ 7 + 2 ^m+^ _|_ q 2 (m-n+fi) _ ]_)) _ q 2 (m-n+fi) 

+ (l-9)2(l-g7+2^ 2 ) f = °' 
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The polynomials {p n {x] q a , g 7 |g)} satisfy the following difference equations 

Pn{x;q a ,q' y \q ) ~ Pn{qx\ q a , q 1 \q) 

g 1 " n (l - q n )(l - q a +'r+ n - 1 ) 


1 - q a+1 

Dq-\ x (w(x- q a , g 7 |g)p n (x; q a , g 7 |g) = 


xp n -i (x; q a+1 , q 1+1 \q ), 

(1 — q a )w(x ; q a ~ l , g 7_1 |g) 


( 6 . 2 ) 


p n+ i(x;<? a \g 7 ^g), 


where 


w(x;q a ,q^\q) = 


(qx-,q)o 0 x c 


(6.3) 


q l ~ n (\ — q n )(l — q m +P+ 7-1' 


Oi^k) 


(g 7+1 x; g)oo' 

Theorem 6.4. For /3,7 > — 1 we have 

M rn,n\ z l. « 2 |?) - MWg\ Zl , qz 2 \q) = - 1 _ qm _ n+ p 

for m > n, and 

D q -I,*! (10(21^2; ? 7 k)^n ) ( 2 ; i ; ^ 1 ?)) 

_ (l-g">-"+^) W ( g i g2 ;^, g 7- 1 | g ) (j3|7 _i) 

^m,n +1 ( 2 1 ) 2 21 <?) > 

z%- m D q -i, Z2 (z™~ n w(ziZ 2 ] g 7 |g)M^)(zr, z 2 |?)) 

- - Z29 m-n+^-l (1 _ ?) - M -,n+l 

/or m > n + 1. 

The above recurrences are obtained by applying the definition of Mm,n\zi, z 2 \q) to (6.2) 
and (6.3). 


7 Polynomial solutions to differential equations 


In this section we study polynomial solutions to partial differential equations. We are looking 
for polynomial solutions to the second-order partial differential equation (3.15). The results of 
this section should be contrasted with those in [9]. 

OO 

Let / = Y1 a j,kZ{ z 2 and substitute in (3.15). 

m,n=0 

Theorem 7.1. The partial differential equation (3.15), namely, 
dd Zl ddz 2 f + ddz ^ = ~ n f 

has a power series solution 


f(zi,z 2 ) = Y aj,kz{z 2 
j,k =0 

if and only if 

OO 

f(zi,z 2 ) = Y a it 


(7.1) 


n\z\ 


3 =1 


{P + 3 + l)i 


-L^ +3 \ziz 2 ) + Y «o ,kZ 22 Fi 


k =0 


—n + k, 1 
k + 1, j3 + 1 


Z\Z 2 
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Proof. Substitute the power series (7.1) for / in (3.15) and equate coefficients of like powers 
of z i and Z 2 to find that 


k — 1 — n 

aj ’ k = w+j) ai - 1 ' k - 1 ’ 

We iterate this and find that 


j, k>0. 


( Tl)kQ>j—k, 0 

k\(/3 + j -k + l)k 
{k- j - n)j(k- j)\ 

a O,k~j 


j > k, 
j < k. 


This proves the theorem. 


Theorem 7.2. In order for the equation 

d Zl d z J + (J^~ z *) d *2 f = x f 

to have a polynomials solution in Z 2 , it is necessary and sufficient that A = — n, n = 0,1,2,, 
in which case the function f will be given as in Theorem 7.1. 


8 A biorthogonal system 


Theorem 8.1. Given three sets A, B, C and a function p from C to B. For each b G B, 
let {5, Bs, dp(x)} and {T, JFt-, dv(x)} be two probability spaces such that {ip a (x]b)f(x,b)} a £A C 
L 2 (S, dp,{x)) and {ip c {y)}ceC C L 2 (T,dv(y)), 


(p ai (x;b)ip a2 (x-,b)\f(x,b)\ dp(x) = ( ai (b)5 ai ,a 


and 


J T ‘ t Pc 1 {y)ip C2 (y)diy{y) = y Cl S cl ,c 2 - 


Then 



® ai ,cx (x, y)$a 2 ,c 2 (x, y)dn(x)u(y) 


Cai (p(ci))r] Cl bai ,a 2 8 C1 ,C 2 i 


where 


$a 1)Cl (x,y) = y?a 1 (a;;p(ci))V , ci(y)/(.T;p(ci)). 

Proof. Observe that for each (a, c) G A x C, & a ,c(x,y) G L 2 (S x T, dp,(x)dv(y)) by Fubini’s 
theorem. Then for 01,02 G A, c\,C 2 G C we have & ai}Cl (x,y)& a2 ,c 2 (x,y) G L(S X T,dp{x)dv(y)) 
by applying Cauchy-Schwartz inequality. By applying Fubini’s theorem we get 


ISxT 


® ai ,ci {X, y)^a 2 ,c 2 (x, y)d/j,(x)i/(y) 


= / VaA x \ p{ci))f{x, p(c{))ip a2 {x\ p{c 2 ))f{x, p{c 2 ))dp,(x) / 'ip Cl (y)ifo 2 {y)di/(y) 

Js JT 

= yc^c llC2 / Ta 1 (x-,p(c 1 ))ip a2 (x-,p(c 1 ))\f(x,p(c 1 ))\ 2 dp(x) = ( ai (p(ci))r? Cl 6 ai , 0 . 2 $C! ,C 2 • 

Js 
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For a, b,c,d G (—1,1), the Askey-Wilson polynomials 

q ~ n , abcdq n ~ l , ae l6 , ae~ ld 


p n (cosd;a,b,c,d\q) = 4 0 3 
satisfy 


ab , ac, ad 


-,Q,Q 


where 


p m (x ; a, b , c, d|g)p n (x; a, b, c, d\q)w{x\ a, b, c, d\q)dx 

_ _ 2-K{abcdq 2n - 1 q) 00 (abcdq n ~ 1 ] q) n 5 mtn _ 

{q n+1 ;q)oo(abq n , acq n , adc/ n , bcq n , bdq n , cdq n ; q)^ ’ 

h(x, l)/i(x, y/q)h(x, -1 )h(x,-y/q) 


w(x:a,b,c,d\q) = ,_, 

/r(x, a)h{x , b)h(x, c)h(x, d)V 1 — x 2 

OO 

h(x, a) = JJ (l — 2axq k + a 2 q 2k ), x = cos 0. 

/c=0 

Theorem 8.2. Let 

^k\ x i y; «i> & 1 , Cl, di; a 2 , b 2 , c 2l d 2 \q) = Pj(x ; ai, ci, dig ,7fc+<5 |g)pfc(y; a 2 , b 2 , c 2 , d 2 |g), 

(7,5)/ , , , , , x Pj(x;ai,6i,ci,dig 7fc+5 |g)p fc (2/;a 2 ,6 2 ,c 2 ,d 2 |g) 

u j,fc Tx,y;ai,6i,ci,di;a 2 ,6 2 ,c 2 ,d 2 |g) =- h(x,d ig ^) -’ 

w(x, y, ai,bi, ci; a 2 , b- 2 , c 2 , d 2 \q) 

_ h{x, l)h(x, y/q)h{x , -l)/r(x, -^/g) hfa l)h(y, y/g)h(y,-l)h{y,-^q) 

h(x, ai)h(x, bi)h(x, ci)Vl - x 2 h(y, a 2 )h(y, b 2 )h{y, c 2 )h(y, d 2 )y/l - y 2 ' 

P^k 3 ’ 7 ’ 5 \ x , y,ai,bi,c 1 ,dy,a 2 ,b 2 ,c 2 ,d 2 \q) 


_ Pj (x; ai, bi, cig afc+/3 , dig 7fc+ ^|g)p fc (y; a 2 , b- 2l c 2 , d 2 |g) 
/i(x,cig aA; +/ 3 ) 

^a^.7,5) ^ y;ai)6l)Cl)dl . a2;62)C2)d2 | g ) 

_ Pj{x\ ai, bi, gq ak+ P, dig 7fc+ ^|g)pfc(y; Q2, b 2 , c 2 , d 2 |g) 
h(x,diq^ k+s ) 


and 


w(x, y; ai,by,a 2 , b 2 , c 2 , d 2 |<?) 

_ L(x, l)/t(x, y/q)h(x, -1 )h(x,-y/q) h(y, l)h(y, yfq)h{y,-l)h{y,-Jq) 

h{x, ai)h(x, bi)Vl - x 2 h(y, a 2 )h(y, b 2 )h(y, c 2 )h(y, d 2 )\Jl - y 2 ’ 

where 


x = cos 0, y = cos a, f3 > 0, 


t/ien 



/-i 


m^ <5) (x, y ; ax, &i, d, di; a 2 ,6 2 , c 2 , d 2 1 g)^ 7 ;,^ (x, y ; ai, &i, ci, di; a 2 , b 2 , c 2 , d 2 1 y) 


x u;(x, y; a 1; ci; a 2 ,6 2 , c 2 , d 2 |y)dxdy 
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_ 47r 2 (q 2 fr 2 c 2 d 2 g 2fc ; gW^^d^^ 1 ; q)k _ 

(q k+ 1 ,a 2 b 2 q k ; q) 00 {a 2 C2q k , a 2 d 2 q k , b 2 c 2 q k , b 2 d 2 q k , c 2 d 2 q k ; g) c 

x _ (aifcicidig 2:,+7fc+ ^; g)oo _ 

X (gJ+i, ai b iq i, aiCiqi +k , a 1 d 1 q^+'r k + s -q)^ 
x (aifeiCi(iig JI+7fc+<? ~ 1 ; q)j5^ m 5 k ,n 

X (6ici^'+ fe , 6 1 di^'+ 7fc +< 5 , cidi^'+ 7fc + 5 ; g)oo ' 


and 


”1 r i 


/-iJ-i 


^(a,/3, 7 ,5) ai,bi,ci,di;a 2 , b 2 , c 2 , d 2 \q) 

x Qm’n ' y,S \ x > y; ai, 6i, Cl, di; a 2 , 62, C2, d 2 |g)w(x, y. ai, 61; a 2 , b 2 , c 2 , d 2 |g)dxdy 
__ 4vr 2 (a 2 6 2 c 2 (j 2 g 2fc ; g) 0 o(a 2 b 2 c 2 d 2 g A: ~ 1 ; g) fc _ 

(q k+ 1 1 a 2 b 2 q k - 1 q)oo(a 2 c 2 q k , a 2 d 2 q k , b 2 c 2 q k , b 2 d 2 q k , c 2 d 2 g fc ; <7)00 

x (aibiCidig^+^^+^^oo 

X (gH 1 , ai&igJ, aiCigi+^+Z 3 , aidigJ+ 7fc + 5 ; g)^ 

(QlfelCid 1 g^'+( Q + 7 ) fc +^- 1 ;g) i d i , m 4,n 

X (6iCigJ+ Qfc +/ 3 ,6 1 d 1 g 7+7fc+l5 , cidig-?+(a+7)M-/3+<5 ; g)^' 

Proof. The details of the calculation are 
/•l /•! 


-1 J-l 


^K/3,7,5) (^ y;ai)6l)Cl ^ i;a2 ,& 2 ,c 2 ,d 2 |g) 

x Qm’n' y,S \ x > y; ai, 61, Cl, di; a 2 , 62, c 2 , d 2 |g)w(x, y; ai, 61; a 2 , b 2 , c 2 , d 2 |g)dxdy 
= J pj (x; ai, 61, cig afe+/3 , dig 7fc+5 |g)p m (x; ai, 61, cig an+/3 , d 1 g rn+<5 |g) 

h(x, l)h(x, y/q)h(x, —l)h(x, — x /g)dx 
h(x, a\)h(x, bi)h(x, c\q ak+ P)h(x, dig 7n+<5 )\/l — x 2 

x [ Pk(y ; «2, b 2 , c 2 , d 2 |g)y n (y; a 2 ,62, c 2 , d 2 |g)tc(y; a 2 ,62, c 2 , d 2 |g)dy 


/-1 


_ 27r(a 2 6 2 c 2 d 2 g 2fc ; g)oo(a2^2C2d 2 g fc g)fc4,n _ 

(q k+ 1 ',q)oo(a 2 b 2 q k , a 2 c 2 q k , a 2 d 2 g fc , 6 2 c 2 g fc , 6 2 d 2 g fc , c 2 d 2 g fc ; g)oo 

x [ pj(x ; ai, 61, cig afc+/? , dig 7fc+<5 |g)p m (x; ai, 61, cig“ fc+/3 , dig 7fc+,5 |g) 


/-1 


h(x, l)h(x, y / g)/i(x, — l)/i(x, — s/q)dx 
h(x, ai)h(x, b\)h(x, ciq ak+ P)h(x, dig 7fc+<5 )Vl — x 2 

__ 47r 2 (a 2 6 2 c 2 d 2 g 2fc ; g) 0 o(a 2 b2C2d 2 g fc ~ 1 ; g) fc _ 

(g fc+1 , a 2 6 2 g fc ; g)oo(a2C2g fc , a 2 d 2 q k , b 2 c 2 q k , b 2 d 2 q k , c 2 d 2 q k ; g) c 

(ai6iCidig 2:;+ ( a+7 ) fc+/3+5 ; g)oo 
X (gJ +1 , ai&igJ, aicigi+“ fc +d, aidigi+T'^ 5 ; g)^ 

(aifeiCidig^ + ( a + 7 ) fc +^ +(5 - 1 ;g) i d i , m 4,n 

X (6icigi+“ fc +d, Mig J ' +7fc+<5 , cidigJ+( Q+7 ) A; +/ 3 + 5 ; g)^' 

The u, v orthogonality can be proved similarly. 

Theorem 8.3. Let 


(x,y;oi,4,ci,di;a 2 ,6 2 ,c 2 ,d 2 |g) = 


Pj( x i a U b\, ci, diq' rk+s \q)pk(y, a 2 , b 2 , c 2 , d 2 |g) 
(dig^+V 0 ; g)oo 
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Then 


■fill ai ’ hl ' Cl ’ dl ’ fl2 ’ b2,C2 ’ d 

x Pm,n(x,y;a 1 ,b 1 ,ci,di;a 2 ,b 2 ,C 2 ,d 2 \q)w(x,y;a 1 ,bi,c 1 ;a 2 ,b 2 ,C 2 ,d 2 \q)dxdy 
_ 47 t 

(q m+ 1 ,a 1 biq m , aiciq m , &icig m , biciq m ] q)^ 

(a\b\C\d\ q ^+^+*; (0161 Cl dl ^+ 7 ^+ 5 - 1 . q)m 


{aid 1 q m +^ n+& , Mig m+7n+<5 , ciciig m +T' ri +' 5 ; g)oo 

_ {a 2 b 2 c 2 d 2 q 2n ] q)oo(a 2 b 2 c 2 d 2 q n ~ 1 ]q) n _ 

X (g n+1 , a 2 b 2 q n , a 2 c 2 q n , a 2 d 2 q n , & 2 c 2 g n , b 2 d 2 q n , c 2 d 2 q n \ q)^ ' 

where x = cos#, y = cos^, a, (3 > 0 ami rc(:c, y; a\, bi, c±; a 2 , b 2 , c 2 , d 2 \q) is the same as in 
Theorem 8.2. 

Proof. Observe that the integral 

rl rl 

J J p j!k (x,y,ai,bi,ci,di;a 2 ,b 2 ,c 2 ,d 2 \q)w(x,y,ai,bi,cr,a 2 ,b 2 ,c 2 ,d 2 \q 


equals 


/: 


x Pmfn{x , y; ai,bi, Cl, dl] a 2 ,b 2 , c 2 , d 2 \q))dxdy 

Pj(x] ai,bi,ci, diq' yk+s \q)p m (x] ai,bi,ci,diq in+5 \q) 

h(x, 1 )h(x, sjq)h(x , —1 )h(x, — v dj)dx 
h{x, a\)h(x, bi)h(x, ci)(d\q^ k+s e l6 , d\q^ n+5 e~ l9 ] q) oo\/l — x 2 

x / Pfc(y; a 2) b 2 , c 2 , d 2 \q)p n {y] a 2 , b 2 , c 2 , d 2 \q)w(y] a 2 ,b 2 , c 2 , d 2 \q)dy 


i -1 


2'KS kjn (a 2 b 2 c 2 d 2 q 2n ] q)oc{a 2 b 2 c 2 d 2 q n 1 \ g), 


(<? n+1 ; q)oc(a 2 b 2 q n , a 2 c 2 q n , a 2 d 2 q n , b 2 c 2 q n , b 2 d 2 q n , c 2 d 2 q n ; g) c 

x J p j (x]ai,bi,ci,diq' yn+5 \q)p m (x]ai,bi,ci,diq' yn+s \q) 

h(x, l)h(x, y/q)h{x, —l)h(x, —y/q)dx 
h(x, ai)h(x, bi)h(x, ci)h(x, diq^ n+s )Vl — x 2 
47 T“ bj,mbk,n 

(q m+1 ,aibiq m , aiciq m , biCiq m , biCiq m ; q)^ 
k {aib l cid l q 2m+ ^ +5 ] q)oo{aibiCidiq rn+in+s ~ l ] q) m 

X (aidiq m +’y n+s , b\diq m+ ^ n+5 , cid 1 q m +'y n+s ] q)^ 

{a 2 b 2 c 2 d 2 q 2n ] q)oo(a 2 b 2 c 2 d 2 q n ~ 1 ] q) n 


X - 


(q n+1 , a 2 b 2 q n , a 2 c 2 q n , a 2 d 2 q n , b 2 c 2 q n , b 2 d 2 q n , c 2 d 2 q n ] q)^ ’ 
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